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Abstract. Correlation functions for matrix ensembles with orthogonal and unitary- 
symplectic rotation symmetry are more complicated to calculate than in the unitary 
case. The supersymmetry method and the orthogonal polynomials are two techniques 
to tackle this task. Recently, we presented a new method to average ratios of 
characteristic polynomials over matrix ensembles invariant under the unitary group. 
Here, we extend this approach to ensembles with orthogonal and unitary-symplcctic 
rotation symmetry. We show that Pfaffian structures can be derived for a wide class 
of orthogonal and unitary-symplectic rotation invariant ensembles in a unifying way. 
This includes also those for which this structure was not known previously, as the real 
Ginibre ensemble and the Gaussian real chiral ensemble with two independent matrices 
as well. 

PACS numbers: 02.30.Px, 05.30. Ch, 05.30.-d, 05.45.Mt 
1. Introduction 

There are many applications of Random Matrix Theory in physics as well as in 
mathematics [U El El El E]- Different matrix ensembles describe universal features of 
eigenvalue statistics in spectra stemming of various physical systems. For example the 
chiral (Laguerre) ensembles with two independent matrices [H El [9J and the Ginibre 
ensembles [TUl [EU [121 EE31 EEH ESI US] describe universal properties of the Dirac operator 
and the Hamilton operator with chemical potential. 

To model physical systems, one has to use Hermitian matrices, if there are no further 
constraints. In the case of time reversal invariance, the matrices to be employed are 
real-symmetric or quaternionic self-dual depending on the behavior of the system under 
space rotations. These three symmetry classes are referred to as unitary, orthogonal 
and unitary-symplectic, respectively. Mean values of characteristic polynomial ratios 
are important quantities to characterize those ensembles. The matrix Green function 
[TT1 ITS] , the replica trick [19] as well as the investigation of the sign problem in 
Quantum Chromodynamics (QCD) [2D] are based on such correlation functions. For 
many ensembles with factorizing probability density it is known that those averages 
have a Pfaffian structure [HI El EH]- The kernels of the Paffians are mean values of 
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one and two characteristic polynomials. Thus, all eigenvalue correlations are completely 
determined by the correlations of the lowest order. 

The method of orthogonal polynomials [211 E21 E3] and the supersymmetry method 
[2U [2SJ E] are successful techniques to derive those Pfaffian structures. Recently, 
we presented a new approach to derive determinantal structures for unitary rotation 
invariant matrix ensembles in a unifying way [26]. Here, we generalize this approach 
to ensembles with orthogonal and unitary-symplectic rotation symmetry. This also 
includes real Ginibre ensembles and Gaussian real chiral ensemble with two independent 
matrices for those the Pfaffian structures were unknown up to now. 

We structure this contribution as follows. In Sec. El we outline our approach. The 
general result is presented in Sec. El We give explicit formulae for the real symmetric 
and the Hermitian self-adjoint matrix ensemble in Sec. HI In the same section, we also 
present two lists of matrix ensembles with orthogonal and unitary-symplectic symmetry 
to which our method can be applied. Concluding remarks are made in Sec. In 



Appendix A[ we explicitly derive some of the equations in Sec. |3j 



2. Outline 

We consider averages of ratios of characteristic polynomials over the Hermitian self-dual 
matrices 

The probability density P is rotation invariant and factorizes in the eigenvalues of H, 
that is, in E = diag (Ei, . . . , E N ) eg) 1 2 . We choose k = diag (ku, . . . , Kki, «42> ■ • ■ , K k2 ) = 
diag (ki, K2) in such a way that the integrals are well defined. The matrix 1 2 n is the 
2A" dimensional unit matrix. Changing to eigenvalue-angle-coordinates yields 

/N k 1 J-, \ 2 

II II nE a ) \ a Zl b2) K(E)d[E] (2.2) 
a=l 6=1 ^ a hl > 

with a normalization constant c. The Vandermonde determinant is defined by 

A^v(^)= II (^-^) = (- 1 )^ 1)/2det K" 1 ]i< a ,^- ( 2 - 3 ) 

l<a<b<N 

Introducing Dirac distributions, we extend the eigenvalue integrals to 2A^ eigenvalue 
integrals and have 

z{k)=c hsm, E j+N ) nn ;_ "" a^ , (2.4) 



J'=l 



where 

5{ Ej — E j +A 

7j 3 ~ E 3+N 



g{E h E j+N ) = P(Ej) K 3 ^' . (2.5) 

hi A hlA 
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In the next step we use the method developed in Ref. [26] based on the idea of Basor 
and Forrester [27|. We extend the product of the characteristic polynomials times the 
Vandermonde determinant by the Cauchy determinant 



Ber (fc/ fc )( K ) 



A fc (Ki)A fc (^ 2 ) 

k 

n (ko,i - K b2 ) 

a,b=l 



(_l)*(fc-l)/2 det 



K a l — «62 



(2.6) 



l<a,6<fc 



Notice that the product 

27V k 

Afc(K 2 ) II H^ Ea - K b 2)A 2N (E) = A k+2N (K 2 , E) 

a=l 6=1 

yields another Vandermonde determinant. Then we have 

Ak{Ki)A k+2N (n 2 ,E) 



(2.7) 



Ber 



(2) 
(k/k+2N) 



[kx; K2,E) = ± 



k k N 

n («oi - ^2) n n («oi - E b 

a,b=l o=l 6=1 



±det 



1 


1 


Kol — K 62 




E b 


K b2 


El 


1 



(2.8) 



times the product of g in the integrand (12. 4p . The notation follows the one in Ref. [26] 
where we have also proven the second equality in Eq. (12. 8p . A Berezinian is a Jacobian 
in superspace. The index (2) at the Berezinian refers to Dyson index /3 = 2. The 
function Ber [^/ fc+2 ^ appears by diagonalizing Hermitian supermatrices a = UsU^ with 
a unitary supermatrix £7 G U (k/k + 2iV), i.e. it is defined by 



f(a)d[a] 



f(UsU^\BeT^ /k+2N) (s)\d[s}dfi(U) +b.t. 



(2.9) 



for an arbitrary superfunction /. The measure dfi(U) is the Haar measure on 
U (k/k + 2N) and "b.t." denotes the Efetov-Wegner boundary terms [281 ESI EQ] which 
may arise in such a transformation. 

As already pointed out in Ref. [22] , this intimate relation to supersymmetry allows 
us to refer to our approach as "supersymmetry without supersymmetry", because we 
never actually map the matrix model onto superspace. Nevertheless, we establish the 
previously unknown connection to superspace. 

Integrating over all energies Ej with j > N in Eq. (12. 4p . we obtain 



Z(k) 



det 



K al — ^62 



K 



a-1 
62 



1 



«al — E b 



Er l 



g(E b ,E) 

Hal - E 



dE 



Jg(E b ,E)E a - 1 dE . 



d[E). (2.10) 
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With help of a modified version of de Bruijn's integral theorem [31], cf. Appendix C.2 
in Ref. [26J, we intergrate over the remaining variables and find the Pfaffian expression 



Z{k) 



Ber (?A)( K ) 



=Pf 






1 

Kbl - 


^a2 




1 


F(« a l 


«6l) 




^62 — ^al 


K b2 


-G„( 


«6l) 


M ab 



(2.11) 



We give a detailed definition of the functions F, G a and M ab in Sec. [31 Here, we 
schematically explain what these functions are. The function F is almost the average 
over two dimensional Hermitian self-dual matrices of two characteristic polynomials in 
the denominator. The functions G a are Cauchy-transforms of the moments of P and M ab 
is the anti-symmetric moment matrix of P generating the skew orthogonal polynomials 
of quaternion type [5J. 

Since the Pfaffian determinant is skew symmetric in the pairs of rows and columns, 
we can construct any linear independent set of polynomials in the last columns and rows 
in Eq. (12. lip . For example, the skew orthogonal polynomials yield a block diagonal 
moment matrix M ab which leads immediately to the well known result expressed in 
terms of skew orthogonal polynomials. Here, we leave the monomials as they are and 
use 



Pf 



A B 
-B T D 



PfDPf [A + BD^B 1 



(2.12) 



for arbitrary matrices A, B and an invertible, anti- symmetric, even dimensional matrix 
D. The matrix A has to be even dimensional and anti- symmetric, too. As M ab is even 
dimensional, we arrive at the final result 

c 



Z(k) 



Ber (?/ fc) ( K ) 



xPf 



2iV 



/J K a2 1 ^mn^b2 1 



m,n=l 



^ 2N 

+J2 G m {n al )M^ n K 

Kb2 — K>al J~ , 



1 .-n-1 
62 



m,n=l 



2N 



^61 — fta2 



^2~ l M^ n G n {K bl ) 



m,n=l 



2N 



F(«ol,«6l) +/ u G m (K al )M m lG n (K b i) 



m,n=l 



Ber (?A)( K ) 



=Pf 



K U (K b2 , K a2 ) 


K X2 


^0,2-, K bl) 


—K 12 (K a x, K b2 ) 


K 22 


, K al, K bl) 



(2.13) 



This is, indeed, the correct result which we found without making use of the Dyson- 
Mehta-Mahoux theorem [32] . Although we can employ an arbitrary choice of polynomial 
set, we obtain the skew orthogonal polynomials generated by M ab . Thus, the skew 
orthogonal polynomials result, not an input. 
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We show in the ensuing sections how Pfaffian structures for a wide class of matrix 
ensembles can be obtained in a unifying way. Our method is applicable not only for 
unitary-symplectic symmetry but also for orthogonal symmetric ensembles. We will 
argue that there is no difference between both symmetries in the derivation. Hence, the 
Pfaffian structure of averages similar to Eq. (12.11) is elementary. 



3. Main result 



We consider the integral 



N 



h(z 2 N+l) Y\ d( z 2j-l, z 2j 



2N+1 k 2 

n no*. 

a=l 6=1 



C 2JV+1 



2N+1 ki 

n n (^bi - z a ) 

a=l 6=1 



-A 



2N+1 



[z)d[z] . 



(3.1) 



We choose the functions h and g and the external variables k = diag (kh, . . . , k^i, K12, 
. . . , Kfc 2 2) in such a way that the integral exists. With the two-dimensional 
Dirac distribution h{z 2 N+i) — 8 (z2N+i) and with the function g(z2j-i,z 2 j) = 
Z2j-iZ2jg(z 2 j-i, z 2 j), we regain another important integral 



Z {2N) (k) 



k 2 2N 

N II II ( K a2 - Z b ) 

Ylg(z2j-i, z 2j )^-^ A 2N (z)d[z\ 

j =i ft n ( K ai - 



(3.2) 



a=l 6=1 

which we calculate in the following. 



7 (2N+1) 
y (fci/te) 



As in Ref. [25], we extend the integrand in Eq. (13. ip by yBer^^^/c) and obtain 

N 

«) = / h(z 2N +l)[[g(z2j-l,Z2j) 



Bel \k\/k 2 +2N+l)( Z ) 



C 2AT+1 



J'=l 



d[z\ 



(3.3) 



where we define z = diag (ki; k 2 , z) . The extension was made in the same way as 
described in the previous section. We then use the determinantal structure (I2.8P of 
the square root Berezinian in the numerator for the integration. In |Appendix A.l| we 



explicitly calculate Eq. (13. 3p for odd d = k 2 — k\ + 2N + 1 > using the sketch of Sec. [2] 
and find 

(-l)WjV!PfM (d) 



(k 1 /k 2 ) v 



X Pf 



BeT ?k\/k 2 )( K ) 



K[f(K a2 ,K b2 ] 



-K{f(K a i, K b2 ) 



l<a,b<k 2 



l<a<fci 
l<b<k 2 



K±%(K bl , K a2 ] 



■K-22 i^alj ^bl, 



l<a<k 2 
l<6<fci 

l<a,6<fei 



(3.4) 
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F(« i,«6i) = - 



(2) 
(2/0) 



— (Kal ~ K bl 



g{z 1 ,z 2 )(z 1 - z 2 ) 



G(d)(^al) 

K(d)(« a2 ) 

•^lf(Ko2, «62) 
•^12 ) («bl, K a2) 




(Kal - Zl)(«al - ^)(«61 - 2i)(«61 - 

C 2 

0(21,22) g(z 1 ,z 2 ) 



d[z] 



«al ~~ ^1 K al — ^2 
-.6-1 -6-1 



d[z] 



h(z) 



Kb<d 



{ Kh a2 l )i 



K 



<fc<d 



(d) (K a2 )M (d 5Kf d) (K 62 ^ 
— + K (d) (K a2 )M^Gf d) ^ / , J ) . 



K&1 — K a2 



) + G(d)(Kal)M f 3Gy)(K M 



{d)^{d)\ 



(3.5) 

(3.6) 

(3.7) 
(3.8) 
(3.9) 
(3.10) 



Here, we use the moment matrix 



M 



(d)= 




w a-l 



g{z 1 ,z 2 )z 1 z 
g(z 1 ,z 2 )z2 



o-l „6-l 
z 2 



- h(z)z a - L dz 



Ka,b<d 



Ka<d 



h(z)z b - l dz 



Kb<d 



(3.11) 



of our probability densities h and g. Let ©m be the permutation group of M elements 
and the function "sign" equals "+1" for even permutations and "—1" for odd ones. We 
fix the sign of the Pfaffian for an arbitrary anti-symmetric 2N x 2N matrix {-D ao } by 

1 N 

Pf [D a b]l<a,b<N = 7^7777 ^2 SigI1 ^ II D ^j)w(2j+1) • (3-12) 



2^iV! 

^e6 2 jv j'=i 

The sums in Eq. (I2.13P for Hermitian self-dual matrices are encoded in the matrix 
products of Eqs. (13.81) . (13. 9p and ( I3.10p . Indeed, we see that the sketch of the approach in 
Sec. |2]can readily be extended from this particular case to the quite general integral ( 13. ip . 

The integral kernels ( 13. 8 p to ( I3.10p have a simple relation to the generating 
function ( 13. ip with other parameters than ki, k 2 and N. We identify them with the 
particular cases {k\ = 0, k 2 = 2), (k x = 1, k 2 = 1) and {k\ = 2, k 2 = 0), 



K[l N+3) (K a2 ,K b2 ) = (-i) N+ K^L~ Kb2 z? N+ ^< 



K 



(2N+1) 
12 1 



K, 



(2AT-1) 
22 1 



«61, K a2) 



\N+1 



iV!PfM (2i v + 3) 
1 



'(0/2) 



MPfM (2JV +i) [Kbl ~ K a2) 



K a 2, K b2) 



7 {2N+l) ( , 



'(1/1) 



1 



AT+1 l^al Kbl 



y (2N+l) ( ^ 
(2/0) \ K al,Kbl, 



iV!PfM (2iV -l) 

The normalization constant is defined by the case k\ = k 2 = 0, 



a 



(2N+1) 



A2N+1) 
'(0/0) 



-l)* +1 iV!PfM 



(2AT+1) . 



(3.13) 
(3.14) 
(3.15) 

(3.16) 
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Hence, we plug these relations into Eq. (13. 4p which leads to the result 

(_l)(k2-*?)/4+fci+ijV! [(-l)^Pf M (d) ] 1_(fcl+fc2)/2 



Z (2JV+1)/ s 



x 



(3.17) 



x Pf 



«&2 



w (cZ-2), 
K a2j^( /2) l K a2 5 ^62 J 

[(d-3)/2]! 
[(d-l)/2]!( K62 - Kal ; 



l<a,fe<fc 2 



l<a<fci 
l<fc<fc 2 



[(d-l)/2]!(K W -K a2 ) 

< x 7 (rf+2)/ s 

[(d + l)/2]! 



l<a<fc 2 
l<6<fci 



l<a,6<fei 



When d is odd, fci + k 2 is even. Thus, the Pfaffians are well defined. Equation (13. IT)) 
implies that the correlations for two characteristic polynomials determine all other 
eigenvalue correlations if the probability density has the factorizing structure as in 
Eq. Q2U). 



For the case that k 2 + k\ is odd, we extend the integral 



Z {2N+1) ( , 



- lim 

KQ 2 — >CO 



Z (ki/k2 + l)( K ) 



K, 



2N+1 
02 



(3.18) 



by an additional parameter kq 2 - This trick is similar to the one in Refs. [331 [26]. Defining 
d = k 2 - ki + 2N + 2 > 0, we apply the result (I3~T7D and find 

(_l)(fca+fti+l)/2jy-j 



z '(fci/fc 2 ) W 



Ber S/ fe2 )(«) 



x 



(3.19) 



x Pf 



7(^-2), 

'(o/i) 1 



^62 J 




l<a<fc 2 



l<a<fci 



[(d- 3)/2]! I 1 ''--^ 

#11 («o2, «62) 



ifi2 (k i, K 62 ) 



l<a,6<fc 2 



l<a<fci 
l<b<k 2 



[{d- l)/2]! I 1 
^2 K «2 



#22 («al, K 6 l) 



l<a<fc 2 
l<fe<fci 



l<a,6<fei 



in the limit K02 — > 00. We notice the appearance of one-point functions. Equation (I3.19P 
is the analog for odd k 2 + k\ to the result ( 13 . 1 T)) which is true for even k 2 + k%. 

The results ( 13 . 1 Tj) and ( 13. 19j) are also true for the integral ( 13. 2p . We simply have 
to choose h as a Dirac distribution. This relation is well known [5] for odd and even 
dimensional ensembles over real symmetric matrices or circular orthogonal matrices. 
Since the probability densities g and h are quite arbitrary this result considerably 
extends the one found by Borodin and Strahov [IB] . 

We are also interested in the case of d = k 2 — k\ + 2N + 1 < 0. Employing the 



sketched derivation in Appendix A. 2 we have the result 
{-l) N N\ 



X 



(3.20) 
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xPf 








^62 — K al 



l<a<fci 
l<b<k 2 



-7 {1) I 






i al / 








^(1/0) 



l<a<fc 2 
l<b<fcl 

l<6<fcl 



l<a<fci 



l<a<fci 
K6<-d 



1 

{F(« a i, «6l)}i< 0)6 <fc 1 



l<a<-d 
l<fe<fci 



For the integral ( I3.2p . we have to omit the column and the row with and have to 

replace d by 2N+k% — k\. The matrix in the Pfaffian (I3.20p is, indeed, even dimensional. 
Thus, the expression is well defined. 

The Pfaffian structure of the sparsely occupied matrix (13 .20j) for d < is a new 
result. The row and the column with Z^j ^ only appears for odd dimensional, real 
symmetric matrices. This factor is the Cauchy-transform of the probability density 
itself. The function F is almost the mean value of the two characteristic polynomials in 
the denominator which has to be calculated, too. However, the N eigenvalue integrals 
are drastically reduced to one or two dimensional integrals. 



4. Applications 

In Sec. I4.1[ we apply the general results to two ensembles over real symmetric matrices 
and Hermitian self-dual matrices. We give an overview of applications for ensembles 
which are rotation invariant under the orthogonal and unitary-symplectic group in 
Sec.i2J 



4-1. Rotation invariant ensembles over real symmetric matrices and Hermitian 
self- dual matrices 

We consider mean values of characteristic polynomials for a rotation invariant 
probability density P over the real symmetric matrices Herm (1, N) or the Hermitian 
self-adjoint matrices Herm (4, N), respectively, 

1] det(H - K j2 t lN ) 

<&(«)= J p ( H )^ (4.i) 

Herm(/3,V) 11 det(H — KjiI^n) 

i=i 

The constant 7 equals one for the real case and two for the quaternionic case. Such 
averages were considered before [2U EU [181 ES]. Here, we apply our method to show 
that the Pfaffian structure arises in a purely algebraic way. As far as we know this is a 
new insight. 

The generating function (14. ip is related to the matrix Green function and thus to the 
fc-point correlation functions. These matrix ensembles describe time-reversal invariant 
systems. 
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For the quaternionic case, the diagonalization of H leads to the identification 

5(x 2 - xi) 



c.f. Eq. fl32D, and 

7 (2N) 



P(x 1 )5(y 1 )5(y 2 



(4)^(^7*2)' 



X\ — X 2 



(4.2) 



C 



A' 



with 



" 2(i-l) 



7T 



J r(2j) ' 



(4.3) 



(4.4) 



Hence, we plug these relations into Eq. (I3.17P for the case c = N + (k 2 — fci)/2 G N and 
find 



Z (iV ' 4) M 

^(fci/fcs)^ 



/-»(4) 7 M) 
°Af ^(0/0) 



a 



Ber 



(2) 



(«) 



(4.5) 



x Pf 



<^b2 — ^a2 



° c (0/2) \ Kj a2)l^b2j 

1 ^(4) 7 (c,4) 

U c-l A (0/0) 



^(l/l)\ K ali K b2) 



/ \ (7(0,4) 

(Kftl - K a2j^( /o) 



(4W C +1,4), 



(«61 - «ol)- 



(2/0) 



^(4) 7 (c,4) 
U 0+1^(0/0) 



(«62 - K al)^ (0 /o) 

The indices a and 6 numerate all variables «; and, thus, the upper left block and the 
lower right block are a k 2 x k 2 matrix and a k% x k\ matrix, respectively. Similarly, one 
finds results for the cases of odd k 2 — ki or negative integer 2N + k 2 — k\ + 1 according 
to the equations (13. 19j) and (I3.20p . 

Let N = 2L + x with x £ {0, 1}- The diagonalization in the real case leads to a 
product of Heavyside distributions Q(Ej +1 — Ej), j G {1, . . . , N — 1}, which is equivalent 
to the ordering of the eigenvalues E\ < E 2 < . . . < E^. Let Zj = Ej + iyj. We split the 
product of Heavyside distributions in two products 

N-l l+x-i L 

J] II 9(V-^)I] e ^-%i)' ( 4 -6) 

j=i j=i j=i 

We absorb the second product of Eq. (14. 6 j) into the probability density and define the 
probability densities 

g( Zl , z 2 ) = g( Zl , z 2 ) = P(E 1 )P(E 2 )6(y 1 )5(y 2 )Q(E 2 - E x ) (4.7) 

and 

h(z) = P(E)6(y) , (4.8) 

according to even and odd N. Due to the integration method over alternate variables 
[36] . we identify 



Z (2L+x) (k) 
z, (fci/fc 2 ) \ K ) 



7 {2L+x,l), 



a 



(1) Z (k 1 /k2) H K J 
2L+x 



(4.9) 
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2L+x 



■IV' 1 IT 

j fi ra/2) 



10 



(4.10) 



Again, we rewrite Eq. (I3.17P with help of relation (14. 9 p and have for even c = + — k\ 



(«) 



°2L+ x Z/ (0/0) 



^JBerj^K) 



(2) 



(4.11) 



x Pf 



(Kfe2 — «a2 



°c+ X -2 Z/ (0/0) 



7(c+X,l), 
'(1/1) 1 



K al, K fc2, 



(«62 -K a i)Z 



(c+X,l) 
(0/0) 



^(1/1)' ( K W> K a2> 



(e+x.i) 

(0/0) 



— Aval 



^(1) 7 (5+X+2,l)/ \ 
U;x " (2/0) (™al> rufelj 



c+X 



c. 



(1) 



.(c+X,l) 



c+x+2^(0/0) 



The indices a and 6 numerate all variables k as in Eq. (14.51) . Similar results are valid 
for the cases of odd k 2 — k\ or negative integer 2N + k 2 — k\ + 1. 
Let d = k 2 — k\ + 7iV > 0. The moment matrices 



M, 



(i) 



(d) 



P(E 1 )P(E 2 )(Er 1 E b 2 - 1 - El E 2 )dEidE 2 

-oo<Ei<E2<oo 

for the real case with even d. 



(4.12) 



Ka,b<d 



M, 



(d) 

/ P(E)E b - 1 dE 

,R 

for the real case with odd d and 



(i) 

(d) 



/ P(E)E a ~ 1 dE 



Ka<d 



Kb<d 



(4.13) 



(4) 



6) / P(E)E a+b - 3 dE 



(4.14) 



l<a,6<d 

for the quaternionic case generate the skew orthogonal polynomials, corresponding to 
the symmetry. Considering the structure of the Berezinian, this shows an intimate 
connection between the method of orthogonal polynomials and the supersymmetry 
method. 

The Pfaffian structures (14. 5 p and (14. lip are well known [18] . However, those Pfaffian 
structures involving the insertion of the relations (14 .3p and (14.91) into Eq. (I3.20p are new 
results. They show that something crucial happens when 2N + k 2 — k\ + 1 is negative. 
The derivation of these structures is purely algebraical. Hence, it is independent of the 
probability density under consideration. Even with help of the supersymmetry method 
one could not reduce the number of integrals in such a substantial way. 

Another new insight of paramount importance is that the structures obtained for 
real symmetric matrices and those for Hermitian self-dual matrices have a common 
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origin. In all the other methods EU [18] both cases were considered separately. Our 
method shows that in both cases, in the real and in the quaternionic one, the underlying 
algebraic structure yielding Pfaffian determinants is the same. 



4-2. A list of other matrix ensembles 

We average ratios of characteristic polynomials similar to the type (14. ip where the 
integration domains are matrix sets different from the symmetric spaces. Those matrix 
sets have to be rotation invariant either under the orthogonal group or under the 
unitary symplectic group. For both symmetries we give a list of ensembles to which 
the integrals (13.11) or ( 13. 2 p are applicable. A real-imaginary part decomposition is 
Zj = Xj + ii/j and in polar coordinates it is Zj = rje lipj . Then, the probability densities in 
Eqs. ( 13 .ip and ( 13. 2 p are equivalent to the probability densities in Eq. ( 14. ip after suitable 
changes of variables. The ensembles with orthogonal symmetry are given in table [T] and 
those with unitary-symplectic symmetry are listed in table |2j 

Our method can be applied to the examples in both tables by the same procedure 



as described in Sec. 14.11 One identifies the probability densities g, g and h in Eqs. ( 13. ip 
and (I3.2p with those obtained from the ensemble under consideration. Then, one can 
use the results (l3T7D . (I3~T9|) and (I3T2C1 . 

All random matrix ensembles in tables [1] and |2] have physical relevance. Real 
symmetric random matrices and Hermitian self-dual matrices model Hamilton operators 
of quantum chaotic systems [15]. Also circular orthogonal and circular unitary- 
symplectical ensembles have applications in quantum chaos. They describe Floquet 
operators in periodically driven systems [15]. Real symmetric and Hermitian self-dual 
chiral ensembles are used to model Dirac operators in QCD [38J. Furthermore one can 
consider all these physical systems in the presence of a chemical potential. This yields for 
real symmetric and Hermitian self dual matrices the corresponding Ginibre ensembles 
and for the chiral case the two matrix models. 

The two-dimensional complex Dirac distribution used in table [1] is defined by 

5 {2) { Zl - z*) = 5{x 2 - x 1 )5{y 2 + y x ) = —S{n - r 2 )%i + (p 2 ) . (4.15) 

We use the short hand notation 

1 ± u 2 , 
V, = ^f, (4.16) 

c.f. Ref. [H]. The functions erfc and K v are the complementary error-function and the 
-fT-Bessel function of order u, respectively. The function / is calculated in Ref. |44J and 
given by 



f(x + w) = 2 / exp 



o 



-2t(x 2 -y 2 )-^ 



K u/2 {2t{x 2 + y 2 ))eifc(2Vt\y\)j . (4.17) 



As in Sec. I4.1[ we notice that the Pfaffian structure arising in all those ensembles 
is fundamental. Particularly, on the formal level, the obvious difference between 



matrix ensemble 


probability density P 
for the matrices 


matrices in the 
characteristic 
polynomials 


probability 
densities g(zi,z 2 ) 
and g(z 1 ,z 2 ) 


probability 
density /i(z) 


rpp 1 c;"\rTTlTTlP'^'^^^ , tti 'a 1~T"i- 

ces [31 m M 


J. 1 II 11 . //t C IN 1 

H = H T = H* 


ii 


P(ri 1 P(rn) x 
x«J(y!)(5(y2)0(»2 - 




circular orthogonal 
ensemble pi] 


P (tr C/ m , m G N) 
17+17 = Ijv and 

r/ T - r7 


U and W 


P(e^ 1 )P(e^ 2 )x 
x5(ri — l)8(r 2 — l)x 

Y(-)| //1 n ] 


P{e vp )S{r - 1) 


real symmetric chiral 
(real Laguerre) 
ensemble 

[2U E3 Ml EH] 


P (tr(AA T ) m ,m e N) 
A is a real N x M 
matrix with 
!/ = M - /V > 


AA T 


P(s 1 )P(s 2 )x 

X(xi£ 2 ) (ly - 1)/2 X 

x5(y!)5(y2)0(^2 - a:i) 


P{x)S{y)x^- 1 ^ 2 


Gaussian real elliptical 

ensemble; for r = 1 
real Ginibre ensemble 

pes ma us um im ee] 

nasal [23] 


exp 1 >tiH T H 


X 


H 


jQ exp [— tx 2 ] x 
ie{i,2} 

x^/erfc(v/2(l + r) % )x 

x[6(y 1 )5(y 2 )e(x 2 - xx)+ 
+2i6M(z 1 -z*)e(y 1 )} 


exp(—rx 2 )5(y) 


x exp -- — - — -tiH 2 


H = H*; 
r > 


Gaussian real chiral 
ensemble [91 EH] 


exp [- tr A T A - tr P T P] 
C = A + /iB 

D = - A T + ^ B T 

A and B are 
real N x M matrices 
with v = M - N > 


CD 


Yl exp [-2r]_ Zj ] x 
je{i,2} 

^ ✓ 1 If /iVo \ ^ ✓ 

X N Vf\ 2r i+ Z 3) x 

x[5{y 1 )5(y 2 )Q{x 2 -xi)+ 

+2zS^(z 1 - z * 2 )e( yi )) 


x v l 2 exp [— 2r]_x] x 
xK v/2 {2ri + x)6{y) 



matrix ensemble 


probability density r 
for the matrices 


matrices m the 
characteristic 
polynomials 


1 1 * 1 * j_ 1 * j_ 

probability density 


|— 1 priTi lfioT) 1 f_ r' 1 nnl in q 4~ri f^c^o. 

1 1C1 1111 l/lCLlla OV'll \ l \ lc\ l lllcXLl itco 

pa us! 


P ftr P m m (= 
# = #t 


11 


P(x )5(v )5(v A X2 ' Xl ^ 
X\ - x 2 


circular unitary-symplectical 

pTICpTTinIp lAI 
CllDClllUlC f \ 


?(tr(7 m ,m6 N) 
rrtrr — -fl 


U and W 


P{e^)8{n - l)x 

x5(r, i) S fa-<Pi) 
sm(v?i - (f2) 


Hermitian self-dual chiral 
(quaternionic Laguerre) 
ensemble [2H E3 ES EH 


P (ti(AA^) m ,me N) 
A is a quaternionic N x M 
matrix with N < M 


AA^ 


P{x l )x^ I - N+1 5{y l )5{y 2 )x 

S(x 2 - xi) 
x— - 

X\ — x 2 


Gaussian quaternionic ellipti- 
cal ensemble; for r = 1 quater- 
nionic Ginibre ensemble [221 E] 


I (t + 1) t 
exp 1 >trH T H 

x exp — - — - — -trH 
H is a quaternionic mat 


X 

rix 


H 


exp [— 2r 2 (sin 2 (px + T cos 2 (fx)\ x 
xri sin(2<^i)5(ri - r 2 )<5(y?i + <p 2 ) 


Gaussian quaternionic chiral 
ensemble [8] 


exp [- tr A* A - tr B^B] 
C = tA + fiB 
D = iA^ + jtBt 
A and P are quaternionic 
N x M matrices 
with v = M - N > 


CD 


K 2v {2i 1+ r 1 )r\ u x 
x exp [2?7_ri cos (p\] x 
xr l sinv9i5(ri - r 2 )5((p 1 + <pa) 



Pfaffian structures 



14 



ensembles with orthogonal symmetry and those with unitary-symplectic symmetry 
becomes immaterial in our derivation. The Pfaffian structure is exclusively due to 
the starting point Eqs. ( 13. ip or ( 13. 2p . Furthermore, we expect that the list of those 
ensembles given here is not complete yet. 

5. Remarks and conclusions 

We extended our method [2B] to integrals of the types (13.11) and (13. 2\ which are averages 
of characteristic polynomial ratios. Those integrals have a Pfaffian structure whose 
kernels are averages over one or two characteristic polynomials. This coincides with 
known results for particular matrix ensembles EU [HI [18] and generalizes those to 
the real Ginibre ensemble and the Gaussian real chiral ensemble. Tables [T] and [2] show 
a wide class of ensembles for which our results are valid. Remarkably, the Pfaffian 
structure arising for ensembles with a real structure as well as with a quaternionic 
structure emerges from the same type of integral. Thus, there is no difference between 
both symmetries from this formal point of view. 

We showed that the Pfaffian structure is a purely algebraic property for factorizing 
probability densities. This is a crucial new insight. It shows that supersymmetric 
structures appearing in the integrand can be found without mapping onto superspace. 
Furthermore, they are the ultimate reason for the existence of the Dyson-Mehta- 
Mahoux integration theorem [32]. Our method also explains why one has found Pfaffian 
structures in so many areas of random matrix theory. As we only perform algebraic 
manipulations, we expect that these structures are even more general. This is indeed 
confirmed by the examples showed in the tables [T] and [2J 

Surprisingly, for the case of a large number of characteristic polynomials in the 
denominator, the kernels reduce to one and two dimensional integrals. These integrals 
are the mean value of one or two characteristic polynomials in the denominator over one 
or two dimensional matrices, respectively. Thus in this case, we have drastically reduced 
the number of integrals, even below the number that would result when mapping onto 
superspace [HU |35j H7]. This new insight shows that there are two regimes depending 
on the number of characteristic polynomials. In the case h% — k\ + 2N + 1 > the 
matrix in the Pfaffian determinant is up to few entries fully occupied whereas for 
^2 ~ k\ + 2N + 1 < we have to take a Pfaffian determinant of a sparsely occupied 
matrix. 

In short, Pfaffian determinants stem in our method from purely algebraic 
manipulations. This is the reason why our results are so general. No integration has to 
be performed. The Pfaffian structures are already contained in the initial integrand. 
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Appendix A. Details of the calculations 

In Appendix A.l we carry out the integrals in Eq. ( 13. 3 p for the case k 2 + 2N + 1 > k\. 



We derive the other case k 2 + 2N + 1 < ki in |Appendix A. 2 
Appendix A.l. The case k 2 + 2N + 1 > k\ 

Let d = k 2 — ki + 2iV + 1 > be odd. We are interested in the integral 

JV 



J h{z 2N+1 ) f[g(z23-U Z2j)\j^ fl 

£2N+1 3 = 1 



/k 2 +2N+l) 



Az)d[z] = (-l) fc 2(fc2-l)/2+7V x ( A _-g 



X 

C 2iV+l 



h{ z 2N+i) Yl g( z 2j-U z 2j) det 
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r i i 




r i i 
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l<b<k 2 




1<&<2JV+1 



1 0/ J Ka<d 1 D J Ka 



l<fe<2iV+l 



l<b<fc 2 

The first step is the integration over all variables Zj with an odd index j. Thus, we have 

N 



I h(z 2N+1 )f[g(z 2j „ h z 2j )^Berfl /2N+1+k2) {z)d[z\ = (_i)M*.-i)/!H-n x (a . 2) 
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Ka<N I C 
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3-1 



Kb<d 
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d[z\ 



K-bl — Z a 

We perform the last integrals with help of a modified de Bruijn's integral theorem [3"Tjl2"o] 
and find 

r N 1 

/ /i(^ + on^(^-i,^ov Ber S/2iv + i +fc2 )( 5 ) rf N = {-i) n+x n\ x 

3=1 

— j {K (d) (Kaa)} 

K bi — K a2 J l<a<fc 2 1 W J l<a5 
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C 2JV+1 



X Pf 
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K al ~ K b2 J l<a<fci 
l<fc<fc 2 
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(A.3) 
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with the matrices defined in Eqs. ( 13. 5113 Finally, we extract the matrix from 
the Pfaffian by inversion, see Eq. (I2.12p . and arrive at Eq. (13 .4p . 

Appendix A. 2. The case k 2 + 2N + 1 < k\ 

Let d = k 2 — ki + 2N + 1 < be an arbitrary integer. Then, we calculate 

N 



C 2N + 1 



^(^iv+i)n^( 2;2 J- i '^)v Ber S/fc 2 +27v+i)( 5 ) rf w 
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Q2N+1 



x det 
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r i i 
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w K a i - Z b ^ 
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(A.4) 



As in Appendix A.l , we integrate first over all variables with an odd index. This yields 

N 

1 



C 2iV+l 



h{z 2N +i) \\g(z 2j -i,z 2j )\J Ber g 



^ kl ( kl -l)/2+ kl - k2 -l J ^ 



K bl — K a2 J l<a<fc2 
l<fe<A;i 



Rr 1 } 




d[z] . 



(A.5) 



Again we use the modified version of de Bruijn's integral theorem and obtain Eq. (I3.20p 
up to the Berezinian in the denominator. 
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